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We propose an extended Nambu-Jona-Lasinio model to include the heavy mesons with the heavy 
quark symmetry. The quark current-current interaction is generalized to include the heavy quark 
currents. In order to comply with the heavy quark spin symmetry at the heavy quark limit, the 
quark mass dependence of the interaction strength is introduced. The light and heavy pseudo-scalar 
and vector meson, their masses and the weak decay constants, are calculated in the unified frame. 
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I. INTRODUCTION 



<n : 

In recent years, some exotic hadron states have been observed in experiments. Many of them cannot be explained 
easily as the conventional quarkonia. A possible interpretation is the hypothesis of molecular state Many of 
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the studies on the exotic states were based on the heavy quark effective theory (HQET). Recently the chiral quark 
^"V model have been used in solving the molecular state [3, Q • Currently, the most difficulty to identify a molecular state 



is the uncertainty of the parameters about the interaction strengths and the form factors. 

On principle these parameters can be calculated from QCD on the quark level. However, in the low-energy region 
where the QCD perturbation method fails we have to rely on effective theories. Among them, the Nambu-Jona- 
Lasinio (NJL) model @, 0] was widely used to investigate many low-energy hadron problems related to the QCD 
symmetries in a simple way [jj-tloj. 
i-C , By means of Dyson-Schwinger equation (DSE) , the dynamic quark mass is generated from the spontaneous chiral 
symmetry breaking. After solving the Bethe-Salpeter equation (BSE), pseudo-scalar mesons are obtained as the 
Qh, Goldstone bosons |g, 0| . Other mesons such as vector mesons and axial- vector mesons were included by introducing 
more chiral invariant interactions [llT - fl3l | . Also the model was extended to comprise the strange flavor [Til UH . A 
bosonization technique was also developed (l6j | and many works were done along this approach [17l - l20l |. 

Because of the QCD color coulomb interaction, a heavy quark spin symmetry is reached in the heavy quark limit 
t-H ' that the dependency of hadronic matrix elements on the orientation of the heavy quark spin vanishes [2l| . From the 
• heavy quark symmetry the HQET formalism was developed (for a review, see Refs. j22l. \25§). 

Some efforts were made on studying the heavy mesons within the NJL model 0, [23- The bosonization technique 
was used in these studies to obtained the meson Lagrangian in HQET. Using the heavy quark propagators in the 
heavy quark limit, the DSE+BSE approach was also been used to calculate heavy meson observables [2q. \2l\. 

In the NJL model study, the color-octet vector current interaction (^A^7 At ^)( , 0Ap7 A 'V') was widely adopted since it 
is closely related to the QCD interaction. In many DSE+BSE calculations such as in Ref. [28j . the interaction between 
two quarks was assumed to be intermediated by the gluon with a complicated effective propagator. So the color-octet 
vector current of the quark should be dominant. The DSE+BSE calculation using the gluon propagator was also 
performed in the heavy meson case [29]. If we naively treat the gluon propagator as a constant in the coordinate 
space, we would obtain an NJL model with the color-octet vector current interaction. In the heavy quark limit where 
the heavy quark mass mq tends to infinity, we will show that the heavy quark spin symmetry is valid only for the 
jj] , color-octet vector current interaction. 

Other contact interactions such as the color-octet axial-vector current interaction ({jj A° 7„ 7s ' ) j^ P ^cl^l^i 3 ) are 
needed to give a more comprehensive description of the light flavor mesons such as the p meson [llTjla ] ■ We will show 
that the heavy quark spin symmetry would not be reached if these interactions exist in the heavy quark limit. To 
maintain this symmetry, these interactions should be considered as higher order terms and should be 1 /toq suppressed. 
This is critical to extend the NJL model to include heavy quark flavors. 

We will extend the NJL model to include the heavy quark flavors. The typical approach of DSE+BSE will be used 
to obtain properties of heavy mesons. In this way, we can calculate the mass splitting between the pseudo-scalar 
mesons D (or B) and the vector mesons D* (or B*) which is the effect of finite heavy quark mass according to the 
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heavy quark expansion. Due to the fact that the heavy quark masses are far beyond the NJL cutoff scale, the usual 
4-dimensional cutoff is not appropriate here. We will use the 3-dimensional cutoff following Refs. |5. ItI lllj. 

In the next section, we will generalize the NJL interaction to include the heavy quark flavor and derive the mass 
dependence of the coupling strength parameters according to the heavy quark spin symmetry. In Section IIII1 we will 
give a brief account of the DSE+BSE formalism to treat the quark and meson states. In Section HVi we will take 
the heavy quark limit and demonstrate the heavy quark spin symmetry. In Section [V[ numerical calculation will be 
performed and the result will be compared to the empirical data. Finally we will give a brief summary. 



II. NJL INTERACTION WITH HEAVY QUARK SYMMETRY 

In many NJL studies, when dealing with the three light flavors q = u,d, s, the interaction was taken to be the color 
current interaction 

£4 = G v {q\ a cl ,q) 2 + G A (q\ a cl , l5 q) 2 . (1) 

The interaction maintains the Ur(3) <Ei Ur(3) <g> SUc(3) symmetry. Here we will not consider the 6-quark interaction 
which was used to deal the ?7a(1) anomaly, since we will not concern the anomaly here and the contribution of the 
anomaly term is small [15| . After a Fierz transformation, we can get a Fierz invariant interaction 

^ = l Gl E [(« A H 2 + (9*75 A'/ g) 2 ] - ^G 2 ]T [{q\)l,qf + (gA} 7 M759) 2 ] + color-octet terms, (2) 

i=0 i=0 



G\ — Gv — Ga, Gi — Gy + Ga- (3) 



where 



Here the Aj's are the flavor Gell-Mann matrices with Ao = y 1 1. The color-octet terms do not contribute to the 
DSE+BSE calculation of the meson. 

In Ref . [24] , where the heavy flavors Q — c,b were considered, only the color-octet vector interaction 
(gAp7 M g)((5Ap7' i Q) was considered. In Section ITVl we will show that only this term in Eq. (JT|) respects the heavy 
quark spin symmetry in the heavy quark limit. 

The color-octet vector interaction is however not enough to describe the light flavor mesons such as the vector 
p meson. We will also show that the heavy quark spin symmetry would not be reached if other than the vector 
interaction exists in the heavy quark limit. To consistently describe the light sector and the heavy sector of the meson 
system, we assume that the NJL interaction is originated in the color-octet vector current. Other currents appear as 
higher order correction in some series expansion and thus should be suppressed by the l/m q factor if the expansion 
is taken with respect to the constituent quark mass m q . According to this thought, we modify the NJL interaction 
Eq. (HJ) to 

U = G V (q\ a c ^q)(q'\ a c ^q') 2 + -^—{q\ a c l^W '^cl^') + J^i&cWvW *cW<l')- (4) 

Here we can take the light and heavy quarks into a unified frame q, q' = u, d, s, c, b. hi and hi are dimensionless 
parameters. m q and m q ' are the constituent masses of the quarks involved in the interaction. 
The Fierz invariant interaction Eq. ^) has the interaction strengths 

Gl=Gy + _!____:, G2 = Gv + h±J}l. (5 ) 

m q m q > m q m q i 

We notice that G\ is tightly related to the quark constituent mass which is the dynamical one generated from the gap 
equation (see Eqs. (0) and (fTUf in Section Hill) . If G\ depends on the constituent quark mass, the gap equation will 
change radically. As shown in FIG. [TJ when hi — hi, the case of usual NJL model where G\ = Gy is independent 
on the constituent quark mass, the gap equation has the two solutions corresponding to two chiral phases which is 
believed to exist in the QCD chiral limit: a Wigner solution at m q = and a chiral symmetry breaking solution at 
m q 7^ when the coupling is large enough. 

However, when h± ^ /12, the gap equation reveals a singularity at m q = 0. Hence, there is no chiral phase (Wigner 
solution). So we must set h\ — hi — h and the NJL interaction turns to be 

A = Gv^Xbi^f + — — [{^X a cl^) 2 + (j\ a cl»l^) 2 ] • (6) 
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FIG. 1: The gap equation where the zero points are the solution of quark mass. The cutoff is taken at A = 750MeV. 

= 0.65, hi = s/12. The dimensionless parameter gv is defined as gv = GvA 2 . We show two typical situations in the figure: 
gv — gc — 9ir 2 /l6Nc where the strength is critically not enough to break the chiral symmetry; and gv = 2.5 with the chiral 
symmetry breaking solution. 



After the Fierz transformation, we obtain the relevant Fierz invariant interaction: for the light sector 

Cf = \Gv [(qX}qf + (qtJ 5 X)q) 2 ] - \ ( G + -^—) [(qX^q) 2 + (qXj^q) 2 ] , (7) 
and for the heavy sector 

£f = \Gv [(Qq)(qQ) + (QH 5 q)(q*l 5 Q)] - \ ( G + -^—) [{Q^q^Q) + (Qj^ 5 q)(qj^ 5 Q)] ■ (8) 

y y y TTiqiriQ j 

Because the difference between the constituent mass and current mass of heavy quark should be small, we will 
ignore the dynamical effect upon the heavy quark mass in the calculation of heavy mesons. 

III. BETHE-SALPETER EQUATION AND MESONS 

Now we will give a brief account of the method of DSE and BSE used in our calculation of meson states. Throughout 
this section we will use Eq. ([2]) as a general form of the NJL interaction. 

The Dyson-Schwinger Equation (DSE) is used to obtain the dynamical quark mass m q . The self-consistent gap 
equation derived from DSE reads 

m q = m° q + E„ (9) 

where m q is the current quark mass and T, q is the quark self energy 

.32 f d 4 p 16Gx 
~ tT -q = *y GlTr J 7^)4 S i(P) = ^- m qh(m q ), (10) 

where S q (p) is the quark propagator. The expression of the integral I\{m q ) is given in the Appendix. 

We use Bethe-Salpeter equation (BSE) to obtain the meson mass and amplitude. The total quark anti-quark 
scattering amplitude is obtained from the ladder approximation. We decompose the amplitude into different Lorentz 
structures [151 ] ■ The relevant amplitudes are 

7p S =T PP (i"f 5 \i ® ij 5 Xj) + T AP (-i$"f 5 \i ® ij 5 Xj) + T PA (ij 5 Xi ® i$jsXj) + T AA {-i^y 5 Xi ® i$jsXj), (11) 
T v ^wW^Xi <8 fu*j), (12) 

where = q^ / ^fq 2 , f]^ = — q tJl q p '. In the ladder approximation, we need only calculate the loop integral 
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which can also be decomposed to 



Then we have 



where 



+ JpA{ilbK ® i|7sAj) + ^iA(-*l75Ai ® i^ 5 Xj), (13) 
^ =»r^(7 M Ai ® 7^). (14) 



16G i/- ■ x , ^ 16Gi., 

ft = — ^— («75Ai ® «75Aj), As =— ^— (A; ® Aj) 

8G 8G (~^) 
A'a = ^-(7^75Ai <E) 7v7bAj), A"v = g - (7aA ® 7f Aj). 

The integrals J^s are defined in ref. [l5| and formulae with 3-dimensional cut-off are collected in the Appendix. 
The meson mass tum is determined by the pole of the amplitude, 

Det(l-JK")| 2 =0. (17) 

To calculate the weak decay constant of a pseudo-scalar meson, the quark-meson vertex is obtained by expanding 
the scattering amplitude near the meson pole. For a pseudo-scalar meson P, which could be 7r, K, D, or B, the qqP 
vertex reads 



V P {p) = i l5 \ l g P (p 2 ) f ~g P (p 2 ) , (18) 

m q + m q > 

where 

& = (t?) K P (l-J AA K A ), (19) 

m q + m q i K A J PA 

9p =— --, j 7T9P, (20) 

mp 1 - JaaK a 



where D — Det(l — JK). The pion decay constant is given by 



"l'/!t)»-' , 7 5 y9(0)h-(p))=i/^%. (21) 



Similar result holds for the kaon decay constant. In the heavy quark case the decay constant is given by 

(0\q(0h^Q(0)\H( P )) = iF H p», (22) 

where H could be D or B. 

IV. HEAVY QUARK LIMIT 

In this section, we will discuss the heavy quark limit. After the Fierz transformation, the relevant interaction 
between a light quark q and a heavy quark Q in a heavy meson is written in the form 

£f ' = |Gi [(Qq)(qQ) + (Qil5q)(qil5Q)] ~ \g 2 [{Ql^iq^Q) + {Q W q){q W Q)] ■ (23) 

Consider the heavy meson at rest, q — muv, v — (1,0,0,0). In the heavy quark limit one assumes that the mass 
difference between ran of the heavy meson and ttiq of the heavy quark is a small quantity Iq, 

m H =m Q +l . (24) 
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The heavy quark momentum p is expanded around the heavy meson momentum q as p — q + k, where k is assumed 
to be far smaller than m,Q. Then the propagator of the heavy quark reduces to 



l 4 + 1 



+ 0- m Q 2(k ■ v + l ) ' 
The expression on the right hand side is independent of rriQ . The BSE loop integrals reduce to 



(25) 



We find 



Jpp = 2iNctr I wr^ 5 $-l q + * e i75 2( k . ttl + ' (26) 

JpA = 2iNct ^! ^^j^tJ-^) 2{k . ttl + u) ' (27) 

/d 4 fc 1 4> + 1 

(2*)«f-m, + i e ils 2<*. „ + („ + >*)• (28) 

J SV = 2iJV c tr„- / || g_^ + i ^ 2(t .^ + ie) . (29) 

^-^/^I^jJ^Tw- (30) 

A l\T f ^ » k ■ V — TYlq . 

J PP = J PA = UNo J j^y 4 {k 2_ m 2 + ie){v . k + lo + i€) > ( 32 ) 

/d k k ■ V + TYlq 

(27T)* (fc»-m» + i6)(t,.Jfe + Zo + i C )' (33) 

J - = W {V-ml) {v .k + l 0) ' (34) 

j -= 4 ^y — (p-^K,. fc+ M — ■ (35) 

After further decompositions J V y = J$ v (g» v - v»v v ) + J^yV^v" and J» AA = J AA (g q ^ - v»v v ) + J 1 AA v^v v , we have 

T r . , T /" d 4 /c k ■ v + m„ 

^- w °/g? »--ix.. t V (36) 

■^-W#(F^5^FW- (37) 

t l . /" d 4 fc k ■ v — m a 



J AA - AiN c I {k 2_ m 2 + ie){v . k + lo+ie) > ( 38 ) 

i d 4 /c — fc • t> — m g 

J ^ = 4 ^ C / J2+Y {k *-ml+ l e){vk + l + l e) ' (39) 



Thus, in the heavy quark limit 



Jpp{h) =Jp A {h) = JjLi(io) = -Jvvih), (40) 

J550o) =Jsv(l<>) = Jyy(M = -JliOo). (41) 

For a pseudo-scalar meson, the mass equation Eq. (I17[) turns to be 

(1 - Jpp((? 2 )Xp)(1 - J L AA (q 2 )K A ) - J 2 pA (q 2 )K P K A = 0, 

which reduces to 

l-(K P + K A )Jpp(lo)=0, (42) 
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TABLE I: Numerical results of the meson masses and decay constants. The cal. I column: results with the NJL interaction 
Eq. ([6|. The cal. II column: results with the interaction Eq. ([T]) for the light meson sector and the interaction Eq. (|23[) for the 
heavy meson sector. The experimental data are taken from Ref. [30( except for Fb and F B which are taken from the lattice 
calculation in Ref. {3lj ] (see also [U H3]). 





cal. I 


cal. II 


exp. 


m«(MeV) 


392 


389 




m s (MeV) 


542 


540 




m^(MeV) 


139 


137 


135/140 


m K (MeV) 


496 


496 


494/498 


U(MeV) 


91.5 


86.5 


93.3 


/if (MeV) 


97.9 


88.6 


114 


m p (MeV) 


771 


775 


775 


m K » (MeV) 


918 


905 


892 


m D (GeV) 


1.87 


1.86 


1.86/1.87 


m Ds (GeV) 


1.95 


1.95 


1.97 


mz,«(GeV) 


1.99 


2.07 


2.01 


m D; (GeV) 


2.12 


2.20 


2.11 


m s (GeV) 


5.28 


5.28 


5.28 


m Bs (GeV) 


5.37 


5.37 


5.37 


m s «(GeV) 


5.31 


5.36 


5.33 


m fl; (GeV) 


5.42 


5.47 


5.42 


Fn(MeV) 


139 


123 


207 


Fds (MeV) 


147 


129 


258 


Fb (MeV) 


96.7 


87.0 


190 (lattice) 


F Ss (MeV) 


107 


91.8 


231 (lattice) 



in the heavy quark limit. The mass equation of the vector partner is 1 — Jy V (q 2 )Kv — 0, which leads to 

l + K v J PP (l )=0, (43) 

in the heavy quark limit. If G2 = G\, the mass equations of the pseudo-scalar meson and the vector meson are 
identical, and the heavy quark spin symmetry is obtained. Otherwise, if Gi 7^ Gi, the mass of the pseudo-scalar 
meson differs from the mass of the vector meson. Similarly, the masses of a scalar meson and its axial-vector partner 
will be degenerate in the heavy quark limit if and only if G2 — G\. 



V. NUMERICAL RESULTS 



In the NJL interaction Eq. ([5]) , the input parameters are the current masses for light quarks and constituent masses 
for heavy quarks, the coupling constants and the 3-dimensional cutoff. We used experimental data of light mesons of 
m T , rriK, m p , f n to determine parameters Tn^, d , m s , Gy, h and A. Then the experimental masses of niD and rag is 
used to determine m c and rrib- The parameters are 

m° u/d =2.79MeV, m° s =72.0MeV, 

m c =1.63GeV, m b =4.94GeV, (44) 
A =0.8GeV, g v =G V A 2 = 2.41, h =0.65. 

The resulted masses and weak decay constants are show in the cal. I column in TABLE fl] We find that the meson 
mass spectra, both the light sector and the heavy sector, are well fitted to the experimental data. One major difficulty 
is that the calculated decay constant decreases with increasing meson mass while the experimental one increases with 
mass. As already shown in Ref. [HI, the theoretical result of fx is smaller than the empirical data. In the case 
of heavy mesons, the theoretical results are smaller than the empirical ones by almost a factor 2. We notice that 
the decay constant increases with the momentum cutoff parameter A. A possible explanation is that the momentum 



m D (GeV) 



m (GeV) 



FIG. 2: The dependence of heavy-light meson masses on tuq. On the left side, with the interaction Eq. l[6|. On the right side, 
with the interaction Eq. (|23p . The light quark is set to u. 



cutoff in heavy sector is larger than in light sector which reflects the fact that the size of a heavy meson is relatively 
small. 

The dependence of heavy meson masses on heavy quark mass are plotted in FIG. [5J We use H to represent the 
heavy pseudo-scalar meson and H* the heavy vector meson. One can see that when the quark mass tends to infinity, 
the mass splitting between H and H* meson vanishes. On the other hand if we use the interaction in Eq. (|23p and 
keep the parameter G 2 unchanged vs the quark masses, i.e. 

G 2 = G V + = 5.41/A 2 , 
ml 

we observe a mass crossing of the H meson with H* meson mass as the heavy quark mass increase. Beyond the 
crossing point, the mass relation is reverted with the H meson above the H*. the mass curve of H will further reach 
the mass threshold and no H bound state exists beyond. So, a naive generalization of the NJL interaction from light 
quark sector to the heavy is inappropriate. 

As a comparison, we also checked with the interaction of a mass independent vector interaction in heavy meson 
sector, i.e. Eq. (|2l?)) with G\ = G2 = g?,k~ 2 . The interaction in light meson sector is Eq. ([T]) with a different set of 
couplings Gi = <?iA~ 2 and G 2 = g 2 A~ 2 . The parameters are 

m° u/d =3.36MeV, m° =81.7MeV, 

m c = 1.68GeV, m b =5.00GeV, A =0.7GeV, (45) 
gi = 2.52, g 2 =5.82 g 3 =2.53. 

The result is shown in the cal. II column in Tab. HI We notice that the mass splitting between the heavy pseudo-scalar 
meson D (or B) and its vector partner D* (or B*) differs from the empirical data by roughly a factor 1.5. For the 
D and D*, it is 210MeV comparing to the empirical data 150MeV, and for the B and B* , it is 80MeV comparing to 
50MeV. To reduce the mass splitting, one may decrease the coupling <?y. But gy can not be too small, otherwise the 
interaction will not be strong enough to bound the D* meson. 



VI. CONCLUSION 



In this work, we have studied light and heavy mesons in a unified frame with the NJL model. We have followed a 
traditional approach of solving the DSE and BSE. We have used a 3 dimensional cutoff to adequately regularize the 
integrals when heavy quarks are involved. 

We have investigated the heavy quark spin symmetry in the heavy quark limit. We find that, in the heavy quark 
limit, the pseudo-scalar meson and its vector partner will have an identical mass equation only if the NJL interaction 
is a color-octet vector interaction which can be recognized as a approximation of a single-gluon exchange interaction. 

Then we propose an extension to the NJL interaction as in Eq. (|23p which introduces the 1 /m q correction to the 
quark current. The mass dependence suppresses the axial- vector current interaction to guarantee that the heavy 
quark spin symmetry still holds in the heavy quark limit. 
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We have performed numerical calculations to the light and heavy pseudo-scalar and vector meson, both for their 
masses and the weak decay constants. The mass spectra fit the experimental data quite well. But the weak decay 
constants always show a large discrepancy to the experiments. A possible explanation is that the momentum cutoff 
in heavy sector is larger than in light sector, which reflects the fact that the size of a heavy meson is relatively small. 
The issue can be studied using some more realistic interaction other than the contact one. 



Appendix A: The Current Condensates in 3D Cutoff 

In BSE, we need to calculate the loop integral 



/j4~. 



r- 



t 



-r'- 



1 



— m' + ie 



where T and V are the interaction vertices. For pseudo-scalar mesons, we have [If 



(Al) 



Jpp = 7,[h{m) + h(m')] + [(to - to') 2 - q 2 ]h{m, to', q 2 ), 



JpA,!i = 9m ( m + m ') 
(m 2 — m' 2 ) 2 



I - 



(to — m'Y 



I 2 {q 2 ,m, m') + g M — (h{m) - h(m')), 



2<7 2 



T L - 
j aa — 



(h - 1$) - (m + m'fh 



(A2) 
(A3) 
(A4) 



For vector mesons, the loop integral is healed by subtracting a certain term Jy V —> Jy V — Jyyiq = 0) + Jy v (q = 0) 



and one can obtain 
1 



vv 



q 2 



2(to 2 + to' 2 )(/ 2 - 1$) - [3(to - to') 



(7 2 -/ 2 )+4(to 2 -to' 2 ) 2 / 2 ° 



The subtracted term tends to zero when one quark mass tends to infinity. The integrations involved are, 

d 4 p I 



I x (to) = 8iN c 



I 2 {m 1 m l ,q z ) = AiN c 



(2tt) 4 (p 2 - to 2 + ie) 
d 4 p 



1 



(2tt) 4 [(p + \q) 2 - to 2 + ie][(p - \q) 2 - to' 2 + ie] 



(A5) 

(A6) 
(A7) 



and define, 7 2 (to,to') = I 2 (m, to',0) and I 2 = d/ 2 /dg 2 | g 2 =0 - After a calculation, one can find, when (to' — to) 2 < 
q 2 < (to + to') 2 , 



h(m) 



Ac 
4tt 2 



4m- 



4TO 2 , n 

1 ^dn 2 



T( > 2 \ NC f 

h(m,m,q) = -^ J 



(v'A 2 +m a + v'A 2 +m' 2 ) 2 ^\ 2 '" 2+ 2 T "' 2 | ( ™ 2 - m ' 2 ) 2 



(m-j-m') 2 



K 2 — q 2 



dn 2 



(A8) 
(A9) 



in which the A 2 is the 3 dimensional cutoff. The same expression can be applied to the case q 2 < (to — m 



i\2 



The integration involved in , I sector is, 



j(l ) = MN C 



d 4 fc 



k • V — TO 



(2tt) 4 (fc 2 - to 2 + ie)(v ■ k + l + ie) 
Assuming v = (1,0), and integrating out ko below the threshold < m q, one can find, 



Vk 2 + TO 2 + Tflq 



k 2 +TO 2 (,/k 2 +TO 2 -/ ) 



(A10) 



(All) 
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Introducing the 3D cutoff, one can get, 

N c f^+O K + 2m„ / 4^" , 



In which, 

k 2 = 4k 2 + 4m 2 (A13) 
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